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ABSTRACT 

We study neutrino emission from the remnant of an inspiraling binary neutron star 
following coalescence. The mass of the merged remnant is likely to exceed the stability 
limit of a cold, rotating neutron star. However, the angular momentum of the remnant 
may also approach or even exceed the Kerr limit, J/M 2 = 1, so that total collapse may 
not be possible unless some angular momentum is dissipated. We find that neutrino 
emission is very inefficient in decreasing the angular momentum of these merged objects 
and may even lead to a small increase in J/M 2 . We illustrate these findings with a post- 
Newtonian, ellipsoidal model calculation. Simple arguments suggest that the remnant 
may form a bar mode instability on a timescale similar to or shorter than the neutrino 
emission timescale, in which case the evolution of the remnant will be dominated by the 
emission of gravitational waves. 
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1. INTRODUCTION 

Considerable theoretical effort is currently directed 
toward understanding the coalescence of binary neu- 
tron stars. The interest stems partly from the promise 
such coalescence holds for generating gravitational 
waves that can be detected by laser interferome- 
ters now under construction, like LIGO, VIRGO and 
GEO. Additional interest is triggered by the idea that 
neutrinos generated during the merger and coales- 
cence might be the source of gamma rays bursts (see, 
e.g., Paczyhski, 1986; Eichlcr et. at, 1989; Janka & 
Ruffert, 1996). Most of the theoretical work has, to 
date, been performed in a Newtonian framework (see, 
e.g., Rasio & Shapiro, 1994; Ruffert, Janka & Schafer 
1996; Zhuge, Centrella & McMillan, 1996). Only a 
few hydrodynamical calculations have taken into ac- 
count post-Newtonian corrections (Oohara & Naka- 
mura, 1995; Shibata, Oohara & Nakamura, 1997), 
while results from calculations in full general relativ- 
ity are very preliminary (Oohara & Nakamura, 1995; 
Wilson & Mathews, 1995; Wilson, Mathews & Mar- 
ronctti, 1996). 

There are a number of issues, however, that cannot 
be addressed even qualitatively in Newtonian gravi- 
tation. Such issues include the collapse of the merged 
remnant to a black hole. Such a collapse is quite 
likely because the mass of the remnant will proba- 
bly exceed the maximum allowed mass for a single 
neutron star. However, the final fate of the collapse 
is completely unknown if the remnant's angular mo- 
mentum parameter J/M 2 exceeds unity. In this case, 
the remnant cannot form a Kerr black hole without 
first losing some of its angular momentum. In prin- 
ciple, angular momentum could be radiated away by 
the emission of gravitational waves, neutrinos, and/or 
electromagnetic (e.g. magnetic dipole) radiation. 

During the coalescence of the two neutron stars 
some fraction of the kinetic energy will be converted 
into thermal energy. Therefore the remnant is likely 
to be hot and emit thermal neutrinos. During the 
rapid plunge of the binary once it reaches the inner- 
most stable circular orbit (ISCO), the angular mo- 
mentum of the system will be conserved approxi- 
mately, so the remnant will be rotating very rapidly. 
In their fully relativistic treatment of corotating, 
polytropic binaries in quasicquilibrium circular or- 
bit, Baumgarte et. al. (1997a, 1997b) found that, 
for a polytropic equation of state with polytropic in- 
dex n = 1, the total J/M 2 of the binary system at 



the ISCO exceeds unity for all binaries in which each 
member has a rest mass <; 0.94M™ ax . Here M™ ax is 
the maximum allowed rest mass of a nonrotating, iso- 
lated, cold neutron star, and Mq, M and J denote the 
total rest mass, energy and angular momentum of the 
binary system. For more realistic, nonsynchronous 
binaries, J/M 2 is likely to be somewhat smaller, but 
may still be appreciable. For many of these bina- 
ries the rest mass of the merged remnant also exceeds 
M™ ax . Immediately following merger, these objects 
may be temporarily stabilized by thermal pressure 
and angular momentum. Only after some of the ther- 
mal energy and angular momentum have been radi- 
ated away can the star later collapse to a black hole. 
The ultimate fate of such an object can only be es- 
tablished by a three-dimensional, relativistic hydro- 
dynamic numerical simulation which employs a hot 
nuclear equation or state and full neutrino radiation 
transport. While the development of such a code is 
underway by several groups, it is still very far from 
complete. 

Our goal in this paper is to isolate the role of 
the neutrino emission in the early evolution follow- 
ing merger and to estimate its effect on the angular 
momentum J and the angular momentum parameter 
J/M 2 . To do so, we consider a simplified scenario for 
the post-merger neutron star remnant and imagine 
that the remnant evolves to an axisymmetric, quasi- 
cquilibrium state on a dynamical (orbital) timescale 
immediately following merger. Since the merger fol- 
lows a rapid plunge from the ISCO, the merged rem- 
nant forms with the same total mass M, the rest mass 
Mq and the angular momentum J of the progenitor 
binary at the ISCO (here we neglect the small loss of 
mass and angular momentum via gravitational waves 
and gas ejection during the dynamical plunge). Ac- 
cordingly, J/M 2 will typically be near unity. We as- 
sume that the merged configuration quickly settles 
into uniformly rotating, dynamically stable quasiequi- 
librium state and focus on its quasi-static evolution. 
We show that the emission of neutrinos is very inef- 
ficient in carrying off J. Depending on the stiffness 
of the equation of state, this emission may lead to 
either an increase or a decrease of J/M 2 . Also, we 
argue that these objects typically develop a bar in- 
stability on a timescale comparable or shorter than 
the neutrino emission timescale. We therefore con- 
clude that while neutrino emission may be important 
for the cooling of the remnant, it is negligible for the 
emission of angular momentum. 
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The paper is organized as follows. In Section 2 
we estimate the rate of angular momentum loss by 
neutrino emission. In Section 3 we incorporate our 
estimate in a post-Newtonian, ellipsoidal model cal- 
culation, that follows the post-merger evolution of 
the remnant. The equations and formalism for this 
approximation are presented in Appendix A. In Sec- 
tion 4 we compare the relevant timescales for angular 
momentum loss by competing processes, and we sum- 
marize our conclusions in Section 5. 

2. ANGULAR MOMENTUM LOSS BY 
NEUTRINO EMISSION 

Consider a neutrino v with four-momentum p a 
emitted from the surface of a uniformly rotating ax- 
isymmetric equilibrium star into an arbitrary direc- 
tion. A local observer with four velocity u a comoving 
with the surface will measure the neutrino's energy to 
be 

8W (e) = -u a p a = -u*pt - u v Pv . (1) 

Here and throughout the paper we adopt geometrized 
units, c = 1 = G. Since the spacetime is stationary 
and axisymmetric, the components pt and p v are con- 
served along the neutrino's path. Therefore, a distant 
observer at rest with respect to the star's center of 
mass will measure the neutrino's energy to be 
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A distant observer will identify the neutrino's angular 
momentum to be 



d,J = —p v . 



(3) 



Combining the last three equations yields the familiar 
result (e.g., Thorne, 1971) 



dM = fldJ - ^-5W {e \ 
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where 51 = u v /it* is the angular velocity of the rota- 
tion. The luminosity as measured by a local, comov- 
ing observer can be defined as 



L 



8W& 



(5) 



where r is the observer's proper time. Equation (J|) 
can then be rewritten 

M =ttj- —L, (6) 



where the dot denotes a derivative with respect to r. 

Equation (^) is fully relativistic. The relation be- 
tween L and J depends on the characteristics of the 
emission of the neutrinos from the star's surface, as 
well as the relativistic structure of the star. Estab- 
lishing this relationship therefore requires detailed 
numerical models of radiating, relativistic, rotating 
stars. Instead, we adopt a simple, first-order Newto- 
nian description to find the approximate relationship 
given by equation (|l4|) below. We assume that the 
neutrinos are emitted isotropically in the rest frame 
of a local stationary observer comoving with the sur- 
face. An observer in a nonrotating static frame will 
therefore find that, on average, each neutrino carries 
off some angular momentum. This angular momen- 
tum can be found by Lorentz transforming from the 
stationary frame into the static frame. To lowest or- 
der in v we find that, on average, each neutrino carries 
off a linear momentum p^ = vE v — QvjE v , where VL 
is the angular velocity of the star, w the distance of 
the surface element from the axis of rotation, and E v 
is the energy of the neutrino in the comoving frame 
(equation (1)). The corresponding angular momen- 
tum is 

P(p = vj 2 QE v . (7) 

The total rate of loss of angular momentum is there- 
fore (see also Kazanas, 1977) 
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where L is the neutrino luminosity and < vj 2 > de- 
notes the average of zu 2 over the surface of the star. 
We can get a reasonable estimate for this average by 
adopting a rotating, compressible ellipsoid model for 
the star (see Lai, Rasio & Shapiro, 1993, hereafter 
LRS). If we denote the principal axes by ai, et2 and 
as (where a\ = ai is measured in the equatorial plane 
and 03 is along the rotation axis), we find 



< w 2 >=-alf(e), 



(9) 



where 
/(e) 

Here 
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is the eccentricity, and 
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is the surface area of the ellipsoid. Note that /(e) ~ 1 
even for large eccentricities. It is convenient to rewrite 
equation (||) in terms of the angular momentum J and 
the moment of inertia /. The latter is I = 2/5 n n Ma\, 
where dimensionless structure constant of or- 

der unity which depends on the star's density profile. 
For a polytropic equation of state 



P = Kp 



l+l/n 



(13) 



where P is the pressure, p the rest mass density, n the 
polytropic index and K the polytropic constant, K n 
can be derived from the Lane-Emden function (see, 
e.g, LRS). Equation (gh can then be written 
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From equation ([l4j), we find that the angular mo- 
mentum J changes by at most (Huet, 1996) 
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assuming that the merged remnant radiates away all 
of its thermal energy [/hot ~ LAt prior to collapse. 
The ratio Uhot/M following merger can be estimated 
by taking the progenitor stars to be cold, with {/hot = 
0, prior to reaching the ISCO. During the subsequent 
plunge and merger, a part of the kinetic energy will 
be converted into thermal energy by contact shocks 
(Rasio & Shapiro, 1994; Ruffert, Janka & Schafer, 
1996). In the extreme case of a head-on collision, a 
very large part of the kinetic energy will be converted 
into heat. The ISCO, however, usually occurs at a 
very small separation of the stars, and therefore we 
expect that only a small part of the kinetic energy 
will be converted into heat. 

To estimate [/hot, we construct post- Newtonian, 
compressible ellipsoid models of hot, uniformly rotat- 
ing stars, generalizing the cold post-Newtonian mod- 
els of Lombardi, Rasio and Shapiro (1997). Matching 
these configurations to binary models at the ISCO 
calculated by Baumgarte et. al. (1997b), we typi- 
cally find Z7hot/[/ C oid ~ 0.2 — 0.3, where [/ co id is the 
total internal energy, excluding the thermal energy 
(i. e. kinetic energy of degenerate nucleons; see Sec- 
tion 3 and Appendix A). Since for typical neutron 



stars Ucoid/M ~ 0.1, we have Uhot/M ~ few per- 
cent following merger. Hence, from equation (|l^), we 
expect that the relative change of the angular mo- 
mentum due to neutrino emission can be at most a 
few percent. We conclude that neutrinos are very 
inefficient in reducing the merged remnant's angular 
momentum. 

We now focus on the change in J/M 2 associated 
with the neutrino luminosity L: 
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Here we have used equations (|^) and ([L4|) as well as 
T rot = J7J/2. The magnitude of J/M 2 therefore in- 
creases if 
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For a dynamically stable rotating ellipsoid with rota- 
tional velocity v < 1 and T rot /\W\ <; 0.27 (see be- 
low), we have T TOt /M <C 1, u* ~ 1 and /(e) ~ 1, so 
that this criterion reduces to 
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K n > 7- 
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(18) 



The critical value of K n = 5/6 corresponds to a poly- 
tropic index of n cr i t = 0.45 (see, e.g., LRS, Table 1). 
Neutrino radiation thus leads to an increase of the an- 
gular momentum parameter J/M 2 for all equations 
of state with n < n clit . Interestingly, the value of 
"crit — 0.45 is very close to typical values that are 
expected in realistic neutron star equations of state. 
Note that we have derived this numerical value in the 
limit of slow rotation. For rapid rotation, as in the 
numerical example of Section 3, we find that the value 
for Ucrit increases in the ellipsoidal approximation. 

For values of n close to n cr jt, the terms on the 
left hand side of equation ( |l7| ) will always be close 
to zero. According to equation jl^), the fractional 
change of J/M 2 will therefore be much smaller than 
the fractional change of J, which we have estimated to 
be in the order of a few percent at most. Nevertheless, 
it is interesting to note that the emission of neutrinos 
can lead to an increase in J/M 2 , and that for nearly 
Newtonian configurations the sign of the change is 
determined essentially by the stiffness of the equation 
of state. 



4 



20 40 60 20 40 60 0.5 1 0.5 1 

t [s] t [s] t [s] t [s] 




Fig. 1. — Time evolution of a rapidly rotating neu- 
tron star in a post-Newtonian ellipsoidal approxima- 
tion (Model 1, see text). We show the luminosity 
L, the angular momentum J/Mq, the principle axes 
a,\ (solid line in the lower left panel) and 03 (dashed 
line), and the angular momentum parameter J/M 2 . 
Although the merged remnant in Model 1 exceeds the 
maximum allowed mass of a cold, nonrotating config- 
urations, it is stabilized against collapse by angular 
momentum. 

3. A NUMERICAL EXAMPLE 

To illustrate these effects, we dynamically evolve 
axisymmetric post-Newtonian, compressible ellipsoid 
models of hot neutron star merger remnants. We em- 
ploy the formalism of Lai, Rasio & Shapiro (1994), 
supplemented by post-Newtonian corrections (Lom- 
bardi, Rasio & Shapiro, 1997) and thermal contribu- 
tions to the internal energy. The thermal corrections 
are taken into account by naively decomposing the 
polytropic constant K in ( |l3|) into a linear sum of a 
cold and a hot contribution, K = K co \d + -Khot- The 
complete formalism and equations can be found in 
Appendix A. 

For initial data we construct equilibrium configu- 
rations as described in Appendix A. 3. These models 
of hot, rotating neutron stars form a three parame- 
ter family, which can be uniquely determined by the 
central density, the eccentricity and the thermal heat 
content (-Khot/-Kcoid)- We can choose these three pa- 
rameters such that the stars' rest mass Mq, total en- 



Fig. 2.— Evolution of Model 2 (see text). The plots 
are labeled as in Figure 1. Although Model 2 is ini- 
tially stabilized against collapse by a combination of 
thermal pressure and angular momentum, a delayed 
collapse is induced as soon as it radiates away a small 
part of its thermal energy. Note that the neutrino 
luminosity L increases during the collapse. 

ergy M and angular momentum J match those of 
binary configurations at the 1SCO as calculated by 
Baumgarte et. al (1997a,b). Here, we summarize re- 
sults for two different progenitor binary models with 
n = 1 and J/M 2 > 1 (see lines 5 and 6 of Table 2 in 
Baumgarte et. al., 1997b, for the binary parameters). 
The two models differ by how close the stellar masses 
are to the maximum allowed mass of an isolated star 
in spherical equilibrium. We pick these two particu- 
lar cases since their evolution turns out to be quali- 
tatively different. For both models, we fix if C old such 
that the rest mass of the remant is Mq = 2 x 1.5M , 
since observed binary neutron stars all have gravita- 
tional masses close to 1.4 Mq, corresponding to rest 
masses ~ 1.5M®. AT C oid is determined by the rela- 
tion K co id = (Mq/ Mo) 2 / 71 , where the nondimensional 
quantity Mq is given in Table 2. 

For Model 1 (line 5 of Table 2 in Baumgarte 
et. al., 1997b), we have -Kcoid = 208 km 2 . In our 
post-Newtonian approximation, this yields a maxi- 
mum allowed rest mass of Mg nax = 1.92M Q (corre- 
sponding to a maximum allowed gravitational mass 
of M max = 1.71M©) for isolated, nonrotating cold 
stars. These values are consistent with those obtained 
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from recent, realistic nuclear equations of state (see, 
e.g., Pandharipande, 1997). The initial gravitational 
energy of Model 1 is M = 2.75Af and the initial 
angular momentum parameter is J/M 2 = 1.05. The 
remnant has an eccentricity of e = 0.935, a thermal 
heat content K^ot /K co u — 0.25 and a% — 28.0 km. 
The thermal energy in this model corresponds to a 
maximum temperature of ~ 14 McV. 

For Model 2 (line 6), we have i^ co id = 173 km 2 , 
corresponding to Af max = 1.75M and Af max = 
1.56M , so that this model's rest mass exceeds the 
maximum allowed rest mass by a larger amount than 
Model 1. It has an initial gravitational mass of 
M = 2.72M and J/M 2 = 1.0, corresponding to 
a remnant of e = 0.94, -Khot/-Kcoid = 0.33, and 
a\ — 23.3 km. Here the maximum temperature is 
~ 28 MeV, which is similar to typical values found by 
Ruffert, Janka & Schafer (1996). 

We show the dynamical evolution of these two 
models in Figures 1 and 2. Even though the rest 
mass of Model 1 greatly exceeds M™ ax , it is dynam- 
ically stabilized, within our approximations, by its 
high angular momentum, even after all the thermal 
energy has been radiated away. By contrast, Model 2 
collapses after it has emitted part of its thermal en- 
ergy. In Model 1, the angular momentum decreases by 
~ 1.5%, while J/M 2 increases by ~ 0.3%. This shows 
that post-Newtonian corrections in an ellipsoid ap- 
proximation, as well as rotation and deviations from 
sphericity, cause the the critical polytropic index n cr it 
to increase to a value greater than one. In Model 2, 
the changes in J and J/M 2 show the same trends but 
are smaller, because the star emits only a small part 
of the thermal energy before it collapses. In a fully 
relativistic treatment, the collapse would not proceed 
directly to a black hole, because the Kerr limit on 
J/M 2 would be exceeded. Binaries with masses closer 
to the maximum nonrotating mass have J/M 2 < 1 
(see Table 2 in in Baumgarte et. at, 1997b). Their 
remnants could collapse directly to Kerr black holes. 

Note also that the neutrino luminosity increases as 
Model 2 undergoes 'delayed' collapse. This is in con- 
trast to spherically symmetric results for delayed col- 
lapse (Baumgarte, Shapiro & Teukolsky, 1996; Baum- 
garte et. at, 1996), in which the luminosity always 
decreases during collapse. This difference can be un- 
derstood from a simple geometrical scaling argument. 
As explained in detail in Appendix A, we assume for 
simplicity that we can write the thermal energy in the 
polytropic form ehot = nifhotP 1+1 ^ n (equation (A4)). 



In the high temperature limit, ehot 
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so that 



T 4 - pl+V" ^ (a 2 a 3 )-( 1+1 /™). Taking the surface 



area of the ellipsoid to be ~ 47ra 2 we find from (A33) 
that in the diffusion approximation the luminosity 
scales as 

L ~ a\- 1/n a- [1+1,n)/2 (19) 
For spherically symmetric collapse, a% — a 3 and 
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For all n > 3/5 (which accomodates 
most realistic equations of state) , spherically symmet- 
ric collapse therefore leads to a decrease of the lumi- 
nosity: the neutrino optical thickness increases faster 
than the thermal energy increases. If, however, the 
star collapses to a pancake, as in Figure 2, a\ will 
remain finite, and L ~ a 3 ( 1+1 / n '/ 2 j n this case the 
increase of the thermal energy overcomes the increase 
of the optical thickness, and the luminosity increases. 
Whenever the surface of a collapsing star approaches 
a newly formed event horizon, the luminosity will al- 
ways be suppressed by the increasing redshift and 
black hole capture. These effects are absent in our 
model calculations. Nevertheless, our results suggest 
that the collapse of a rapidly rotating, hot neutron 
star following merger may leave a very distinct signa- 
ture in the neutrino signal. 

4. ESTIMATE OF TIMESCALES 

So far we have focussed on the emission of neu- 
trinos and their role in the evolution of the merged 
remnant. To estimate if this evolution is indeed gov- 
erned by neutrinos, we now compare the timescale 
for the neutrino emission with those associated with 
eletromagnetic and gravitational radiation. 

The timescale for the emission of the neutrinos is 
given by the diffusion time, r v . Taking n to be Rosse- 
land mean opacity due to scattering off nondegenerate 
neutrons and protons (see, e.g., Baumgarte, Shapiro 
& Teukolsky, 1996), we find 
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Assuming that electromagnetic emission is domi- 
nated by a magnetic dipole radiation, its timescale 
can be estimated from (see, e.g., Shapiro & Teukol- 
sky, 1983) 
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Clearly, tem is by many orders of magnitudes larger 
than r„. Electromagnetic radiation is therefore negli- 
gible in the early evolution of hot neutron stars, unless 
the magnetic fields are extreme {B ;> 10 15 G). 

An axisymmetric, stationary star does not emit 
gravitational waves. However, if the ratio of the 
kinetic and potential energy t = T Iot /\W\ is large 
enough, the star will develop a bar instability and 
will then emit gravitational radiation. More specifi- 
cally, the star is secularly unstable to the formation of 
a bar mode if t > t scc ~ 0.14, and dynamically unsta- 
ble if t > idyn ~ 0.27 (Chandrasekhar, 1969; Shapiro 
& Teukolsky, 1983; see Stergioulas & Friedman, 1997, 
and Friedman & Morsink, 1997, for a recent relativis- 
tic treatmen t). W ritin g t in terms of the ellipsoidal 



expressions (A10) and (A18), we find 
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Since relativistic configurations typically have M/R ;> 
0.1, this suggests that all neutron stars, rotating with 
J/M 2 ;> 1, have t ;> 0.2 and are at least secularly 
unstable to the formation of a bar mode. This is 
in agreement with the general relativistic, numerical 
models of Cook, Shapiro & Teukolsky (1992). For 
tscc < t < idyn j the timescale for the formation of 
a bar mode driven by gravitational radiation can be 
approximated by (see, e.g., Friedman & Schutz, 1975; 
Lai & Shapiro, 1995) 



T ba 



Note that shear viscosity can also drive a bar mode 
instability, but it is inefficient in hot neutron stars 
with T ~ 10 MeV as is the case here (see Bonaz- 
zola, Frieben & Gourgoulhon, 1996, and references 
therein). Once a bar has fully developed, the grav- 
itational radiation dissipation timescale can be esti- 
mated from the quadrupolc emission formula 
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For t < idyn) the total gravitational radiation time- 
scale tgw ~ T bar + T gra v is therefore dominated by the 
initial growth time Tbar- Since Tbar strongly depends 
on t, this timescale is quite uncertain. However, the 
estimate ( p2]) suggests that for compact stars with 
J/M 2 ;> 1, we have tqw ^ 3 s. We then see that 
the evolution of these merged stars is characterized 
by the formation of a bar mode (see Lai & Shapiro, 
1995) and emission of gravitational waves, which may 
then carry off angular momentum very efficiently. In 
some cases, t > idyn> and a bar may develop on a 
dynamical timescale (~ ft^ 1 ~ few ms; see Rasio & 
Shapiro (1994) for a numerical demonstration). 

5. CONCLUSIONS 

We conclude that neutrinos are very inefficient 
in carrying off angular momentum from hot, mas- 
sive remnants of neutron star binary mergers, even 
if these are rapidly rotating and have an angular 
momentum parameter J/M 2 ;> 1. Curiously, neu- 
trino emission may even increase J/M 2 by a small 
amount. We find that the timescale for the formation 
of a bar mode may, for such objects, be much shorter 
than the neutrino emission timescale. For pure grav- 



(M/R) 2 (m) 6 {M/R) 2 {my 



itational radiation the third term in equations (16) 
and ( |l7| ) vanishes, so that J/M 2 always decreases. 
This suggests that gravitational waves will be faster 
and more efficient in reducing J/M 2 than neutrinos. 
Together, these arguments may have important con- 
sequences for gravitational wave detectors, since they 
suggest that the evolution of mergered binary rem- 
nants will be dominated by the emission of gravita- 
tional waves. In particular, should J/M 2 be larger 
than unity upon merger, gravitational waves seem to 
be the only means by which J/M 2 can be efficiently 
reduced to a value smaller than unity, so that the 
entire merged configuration can collapse to a black 
hole. Otherwise, hydrodynamic stresses will have to 
support or expel some of the matter to allow the in- 
terior regions to collapse. 

These arguments imply that, except for cooling 
and deleptonizing, neutrino emission plays a very mi- 
nor role in determining the mass and spin of the final 
configuration (black hole or rotating neutron star). 
However, because of the role in inducing delayed col- 
lapse in some cases, the neutrinos will have to be in- 
cluded in calculations of the late stages of coalescence. 
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Gravitational radiation will play a very important role 
in determining the final mass and spin of the black 
hole or neutron star. This result strengthens the con- 
clusion that a fully general relativistic description of 
binary neutron star mergers is essential for full un- 
derstanding. 

Note that our arguments apply to massive, hot, ro- 
tating neutron stars in general, and not only to the 
remnants of binary neutron stars. For example, they 
equally apply to newly formed, hot neutron stars in 
supernovae. It has recently been suggested that these 
stars may be stable initially, but may later delep- 
tonizc, undergo a phase transition and collapse to a 
black hole (see, e.g., Brown & Bethe, 1994). Alterna- 
tively, such a delayed collapse might result from the 
loss of angular momentum from a nascent, rapidly 
rotating neutron star which is stabilized initially by 
its high spin. Neutrino emission could, in principle, 
reduce this angular momentum and hence induce col- 
lapse. While for type II supernovae Uhot/M is slightly 
larger than in the examples presented in Section 3, 
our arguments show that the neutrino emission is still 
very inefficient in carrying off angular momentum. 
This fact makes this scenario very unlikely (Huet, 
1996). Similarly, our analysis suggests that neutrino 
emission plays a negligible role in determining the 
quasi-steady spin rate of accreting neutron stars at 
the center of Thorne-Zytkow objects, which have re- 
cently been suggested to be likely sources of gravita- 
tional waves for detection by GEO (Schutz, 1997). 
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A. POST-NEWTONIAN, ELLIPSOIDAL 
MODELS OF HOT, ROTATING NEU- 
TRON STARS 

A.l. The Euler-Langrange Equations 

The total energy of an axisymmetric, rotating, hot 
neutron star can, in a post-Newtonian, ellipsoidal ap- 
proximation, be written as the sum 

M = U + U PN + W + W PN + T, (Al) 



where U is the internal energy, W the potential en- 
ergy, T the kinetic energy, and Upn and Wpn are 
post-Newtonian corrections to the internal and po- 
tential energy. 

Including the post-Newtonian terms not only im- 
proves the solution for relativistic stars quantita- 
tively, but also, more importantly, it changes the 
solution qualitatively. Without the post-Newtonian 
corrections, a graph of the rest mass Mo as a func- 
tion of the central density p c exhibits no turning 
point. By including the post-Newtonian terms, we 
find turning points reasonably close to the turning 
points of the corresponding fully relativistic (Tolman- 
Oppenhcimcr-Volkov) solution (see Lombardi, Rasio 
& Shapiro, 1997). Since these turning points give 
the maximum allowed mass configuration, the post- 
Newtonian corrections allow us to approximate the 
onset of radial instability and to mimic the delayed 
collapse of hot, rotating neutron stars to black holes. 

We assume a polytropic equation of state 

P = Kp 1+1 ' n , (A2) 

and crudely account for thermal pressure by taking 
the constant K to have a cold contribution due to 
degenerate nucleons, K co i d , and a hot contribution 
due to thermal heating, If hot: 

K = K cold + K hot . (A3) 

The internal thermal energy density ehot is therefore 
taken to be of the form 

£hot = nK hotP 1+1 / n . (A4) 

While K co \ d is assumed to be constant in both space 
and time, we allow Khot to depend on time, keeping 
it constant throughout the star. Integrating the en- 
ergy density over the star allows us to write the total 
internal energy in terms of the central density p c as a 
sum of a cold and a hot contribution 

U = U cold + C/ hot - h {K cold + K hot )pll n M (A5) 

where k\ is a dimensionless structure constant that 
depends on the density profile, i.e. the polytropic in- 
dex n. Numerical values for fci can be found, for 
example, in LRS (see Table 1). In terms of the mean 
radius 

R = ( a 2 a3 ) 1/3 , (A6) 
the central density can be written 



8 



where we have defined 



and a contribution due to spin 



6 



47^(^)1' 



(A8) 



and where and £ are the usual Lanc-Emdcn param- 
eters for a polytrope (see, e.g., Chandrasekhar, f 939). 
The internal energy can therefore be rewritten 



l+l/r 



R 3/n ■ 

The potential energy W is given by 



(A9) 



W 



^M|sin^e (1 _ e2)1/6 (A1Q) 



5 — n R e 
3 



5 — n 



/3- 1 / 3 M 5/3 pl/3£^ (1 _ e 2 ) l/6 



Introducing the dimensionless coefficients Ai (Chan- 
drasekhar, 1969, §17) 



li = (1 _ e2)1/2 £^_iz^ 



Al = (l-^)V 2 — — (All) 

and 

A 3 = 2- 2Ai (A12) 
(for ai > 03), the potential energy can be rewritten 

3 M 2 J 



W 



5-n R 2i? 2 ' 



(A13) 



where 

J = 2A 1 a 2 1 +A 3 a 2 3 . (AH) 

The last term J / (2R 2 ) takes into account corrections 
to the potential energy due to deviations from spheric- 
ity. Note that partial derivatives of J take the form 



dj_ _ 1_ 

dcij a, 



(J - a 2 Ai), 



so that 



dW 3 M 2 



CLi Ai 



(A15) 
(A16) 



dat 5-n2R 3 
(see Chandrasekhar, 1969). 

Assuming the star to be uniformly rotating, the ki- 
netic energy T can written as a sum of a contribution 
due to expansion 

T cxp = ^K n Ah{2a\ + al) (A17) 



T rot = \m 2 = J 2 M -5/3 p 2/3 (1 _ ^1/3^ 

2 4k„ 

(A18) 

Here the moment of inertia is given by 
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(A19) 



and structure constants of order unity that de- 

pends only on the polytropic index n (see, e.g., LRS) 
Following Lombardi, Rasio & Shapiro (1997), the 
post-Newtonian corrections can be written 



UpN 



■ M, 



2+1 /n 



^^ 1/n+1/3 ii + 3Ar ( A2 °) 



and 



W PN = -hp 2 J 3 M 7 /3 = -hP 2/3 ^-. (A21) 

Numerical values for the structure constants l\ and I2 
are calculated in Lombardi, Rasio & Shapiro (1997). 
Both post-Newtonian terms are correct in this form 
only for spherical configurations. However, we shall 
assume the slow rotation limit (T ro t/|W| -C 1), for 
which corrections arising from deviations from spheric- 
ity are higher order and will be neglected. 

From these energy contributions we can now con- 
struct the Lagrangian 



L(q i ;q i )=T-U-U J 



PN 



w-w. 



PN, 



(A22) 



where the % are the independent variables and 4> 
(with (f> = fi). The dynamics of the system is then 
governed by the Euler-Lagrange equations 



d_dL_ <9L _ (j 

dr den dqt 



The angular momentum can be defined as 



J : 



dL 

on' 



(A23) 



(A24) 



Since <f> is an ignorable coordinate of the Lagrangi- 



an (A22), the Euler-Lagrange equations imply that J 
is conserved on a dynamical timescale. However, we 
want to allow for the radiation of angular momentum 
on a secular timescale by the emission of neutrinos, 
so we will therefore allow for an external dissipation 
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term on the right hand side of equation (A23) when 

~{J)v (A25) 



dJ 



Note also that when we allow Kv, ot to decrease with 



time as the star cools (see equation ( A39 ) below), the 
system will no longer conserve energy. 

We now find that the Euler-Langrange equations 
yield 



Kf3 l/n^l M Q 



1/n 



i? 3 /"ai 



(A26) 



K 



n 3 / 



-1/3 



5h M 



l/n+l 



10 „o/q h 



M 2 



«3 



— /3 2/3 - 
3 k„ i? 2 ai 



i? 3 /"+ 1 a 1 
3 M 



2n n 5 — n R 3 



K 



R 3/n 



(A27) 



5Zi Mq 



l/n+l 



10 o 2/3 > M 2 



K „ i?3/«+l a3 

5 3 M 



and 



K n R 2 a 3 2k„ 5 



1 ai 



R 3 



a3 A 3 



(A28) 



In the absence of neutrino emission ((■/)„ = 0), these 
ordinary differential equations completely describe 
the adiabatic evolution of a uniformly rotating, ax- 
isymmetric neutron star in a post-Newtonian ellip- 
soid approximation. They represent in this limit the 
ellipsoidal analogue of the post-Newtonian equations 
of hydrodynamics for an adiabadic gas. In order to 
allow for cooling and loss of angular momentum, we 
have to supplement this system of equations with ex- 
pressions for Khot and {J)u- 

A. 2. Neutrino Cooling and Angular Momen- 
tum Loss 

We approximate the thermal energy contribution 
from hot nucleons by 



£nuci — 2 



(A29) 



(strictly true in the limit of high temperatures) , and 
add to it the thermal radiation of photons, electron- 
positron pairs and N v — 3 flavors of (nondegenerate) 



neutrinos 

£rad = (3 + 7N„/8)aT 4 . (A30) 
The total thermal energy density is then given by 



^hot — ^rad 



(A31) 



This expression can be used to estimate the tempera- 
ture T. Adopting the Rosseland mean opacity arising 
from scattering off nondegenerate neutrons and pro- 
tons, the neutrino opacity takes the form 



7(2^) 2 <j SkTY 
20 4m B \m e J 



(A32) 



where me is the baryon mass, m e the electron mass, 
and (To = 1-76 x 10~ 44 cm 2 . The neutrino luminosity 
can then be estimated via the diffusion formula 



L = 



A 



IVe 



racl | 



A (A33) 



where p a = 3Mo/(4ttR 3 ) is the average density. Since 
a 3 < oi, we expect that the "preferred" direction of 
diffusion is along the rotation axis, and therefore esti- 
mate the gradient of e ra d by dividing e ra d by 0,3 . Using 
equation (j^) , the angular momentum loss can now be 
written 

(j) v = < m 1 > Q.L. (A34) 

Next we have to relate the luminosity L to i^hot- 
From equation (^|) we have 



L 



SW (e 
dr 



<T Tt>> 



(A35) 



where the brackets denote an average over the star 
and S is the total entropy of the star. From the first 
law of thermodynamics we have 



Tds 



Pd 



(A36) 



Here e = nKp 1+1 ^ n is the total internal energy den- 
sity, and s is the specific entropy (per rest mass). For 
adiabatic changes, K is constant and Tds — 0. Here, 
K = -Kcoid + ^hot is not constant, however, and we 
find 

Tds = np 1/n dK. (A37) 

In a Newtonian treatment, integration over the whole 
star now yields 



<TdS> 



(Tds)dm = dK np 1/n dm 



d In K J edm = Ud In K. (A38) 
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In keeping with our earlier assumption, we take dK 
to be constrant throughout the star. To incorporate 
the post-Newtonian correction consistently, we sim- 
ply replace U with U + Upn- Inserting this result 
into ( A35| ) then yields 



if hot = — - 



K 



U + Up N 



-L, 



(A39) 



Equa tions ( |A26| ) to ( |A28| ) together with flA34| ) and 
(A39) now approximate the evolution of a cooling, 
rotating neutron star. At each timestep, the luminos- 
ity L can be estimated from (A33). 



Alternatively, we could have derived equation ( A38) 
from observing that, for quasi-static changes of con- 
stant rest mass, the mass-energy of an equilibrium 
configuration changes according to 



dM = fid J + AdC+ < TdS >, 



(A40) 



where we have set u ~ 1, and where for completeness 
we have allowed a vorticity term AdC. Here C is the 
circulation and A is the vorticity of the configuration. 
Following Appendix D in LRS, this change can also 
be written 



dM = 



dM 



dJ 



dJ- 



cci,C 



dM 
~dC 



dC+ 



dM 



dai 



dcti. 
(A41) 

Here the ai are the independent parameters of the 
configuration, which we can take to be p c , e and if hot- 
Note that 

(W =0 (A42) 

Qi,C 



dJ 



and 



dM 
~dC 



= A 



(A43) 



on, J 



(see LRS). Furthermore, since we are considering 
quasi-static changes between equilibrium configura- 
tions, first-order variations with respect to the dy- 
namical parameters vanish, so 



dM 
dp c 



e,K hot ,J,C 



dM 



= 



Pc,K hot ,J,C 



(see Section A. 3). Finally, from (Al), 



dM 



dK ho 



U + Upn 
if 



(A44) 



(A45) 



p c ,e,J,C 

Putting terms together, we therefore find 

dM = CldJ + AdC + (U + U PN )dhxK. (A46) 



Comparing this expression with ( A40| ) shows that we 
can identify 



< TdS >=(U + U PN )d\riK, 



(A47) 



which agrees with (A38) (with U replaced by U + 
U PN ). 

A. 3. Equilibrium Configurations 

A configuration in equilibrium has T exp = 0, and 
the equilibrium relations can then be found, for ex- 
ample, by setting the derivatives of the energy func- 
tional (Al) with respect to the central density p c and 
the eccentricity e to zero: 



= 0. 



dM _ dM 
dp c de 

The first condition yields the virial relation 



(A48) 



l u+ ^±^ UpN + \ w + lw PN + \r n 

n An 6 6 6 



0. (A49) 



Since only W and T rot depend on the eccentricity 
e, the second condition yields in this approximation 
the same result as a purely Newtonian treatment (see 
equation (3.21) in LRS) 



t = 



T 

W\ 



3 



1 - 



e(l-e 2 ) 1 /2 



sm 



- 1. 



(A50) 



Note that t only depends on the eccentricity, and is in- 
dependent of, for example, the mass or the polytropic 
index n. 

For a given central density p c , eccentricity e and 
thermal contribution if hot, equilibrium conditions 
can now be constructed by first finding t from equa- 



tion ( A50). This value can then be inserted into ( A49) 
to eliminate T rot : 

i U+ ^±Jl UpN + lw(l-2t) + ^W PN = 0. (A51) 
n An 6 A 



Writing the individual energies in terms of the rest 
mass Mq and the central density p c , we obtain a 
quadratic equation for M ( 



2/3 



AM, 



4/3 



BM, 



2/3 



, r 2/3 B ( B 2 C 



C = 

1/2 



(A52) 



(A53) 
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where 
A = 

B = 
C = 



hp- 1 " 



1 /3 -l/3 p -2/3™_f (1 _ e 2 ) l/6 (1 _ 2t) + 



5 — n 
3 + n 
3n 



l/n-2/3 



(A54) 



1 



1/n-l 



In ( A53 ), the smallest positive solution is the physi- 
cally relevant solution. In the Newtonian limit, A = 
0, and M is given by (C/B) 3 / 2 (see equation (3.22) 
in LRS). Once Mo has been determined as a function 
of p c , e and if hot, the individual energy terms as well 
as the total energy can be calculated. 

REFERENCES 

Baumgarte, T. W., Janka, H.-T., Keil, W., Shapiro, 
S. L., & Teukolsky, S. A, 1996 ApJ, 468, 823 

Baumgarte, T. W., Shapiro, S. L., & Teukolsky, S. A, 
1996, ApJ, 458, 680 

Baumgarte, T. W., Cook, G. B., Scheel, M. A., 
Shapiro, S. L., & Teukolsky, S. A, 1997a, Phys. 
Rev. Lett., 79, 1182 

-, 1997b, submitted 

Bonazzola, S., Friegen, J., & Gourgoulhon, E., 1996, 
ApJ, 460, 379 

Brown, G. E., & Bethe, H. A., 1994, ApJ, 423, 659 

Chandrasekhar, S., 1939, An Introduction to the 
Study of Stellar Structure (University of Chicago 
Press, Chicago, Illinois) 

Chandrasekhar, S., 1969, Ellipsoidal Figures of Equi- 
librium (Yale University Press, New Haven, Con- 
necticut) 

Cook, G. B., Shapiro, S. L., & Teukolsky, S. A, 1992, 
ApJ, 398, 203 

Eichler, D., Livio, M., Piran, T., & Schramm, D. N., 
1989, Nature, 340, 126 

Fri edman, J. L., fc Morsink, S. M., 1997, submitted, 
( |gr-qc/9706073|) 

Friedman, J. L., & Schutz, B. F., 1975, ApJ, 199, 
L157; Erratum, 1978, ApJ, 221, L99 

Huet, P., 1996, private communication 

Janka, H.-T., & Ruffert, M., 1996, A&A, 307, L33 



Kazanas, D., 1977, Nature, 267, 501 

Lai, R., Rasio, F. A., & Shapiro, S. L., 1993, ApJS, 
88, 205 (LRS) 

Lai, R., Rasio, F. A., & Shapiro, S. L., 1994, ApJ, 
437, 742 

Lai, R., & Shapiro, S. L., 1995, ApJ, 442, 259 

Lombardi, J. C, Jr., Rasio, F. A., & Shapiro, S. L., 
1997, Phys. Rev. D, 56, 3416 

Oohara, K.-L, & Nakamura, T., 1995, Lecture deliv- 
ered at Les Houches School "Astrophysical Sources 
of Gravitational Radiation" (Les Houches, France, 
1995), to be published in the Proceedings (eds. J.- 
A. Marck & J. -P. Lasota) 

Paczyhski, B., 1986, ApJ, 308, L43 

Pandharipande, V. R., 1997, to appear in "Proceed- 
ings of the 10th International Conference on Re- 
cent Progress in Many-Body Theories" , Sydney, 
Australia, ed. D. Nelson, (World Scientific, Singa- 
pore) 

Rasio, F. A., & Shapiro, S. L., 1994, ApJ, 432, 242 

Ruffert, M., Janka, H.-T., & Schafer, G. 1996, A&A, 
311, 532 

Schutz, B. F., 1997, to appear in "Mathematics of 
Gravitation. II. Gravitational Wave Detection", 
ed. A. Krolak (Banach Center, Warsaw) 

Shapiro, S. L., & Teukolsky, S. A., 1983, Black Holes, 
White Dwarfs, and Neutron Stars (New York: Wi- 

ley) 

Shibata, M., Oohara, K.-L, & Nakamura, T., 1997, 
Prog. Theor. Phys, in press 

Ste rgioulas, N., fe Friedman, J. L., 1997, submitted 
( kr-qc/9705056|) 

Thorne, K. S., 1971, in Y. B. Zel'dovich and I. D. 
Novikov, Relativistic Astrophysics, Vol. 1 (Uni- 
versity of Chicago Press, Chicago, Illinois), Sec- 
tion 10.7 

Wilson, J. R., & Mathews, G. J., 1995, Phys. Rev. 
Lett., 75, 4161 

Wilson, J. R., Mathews, G. J., & Marronetti, P., 1996, 
Phys. Rev. D, 54, 1317 

Zhuge, X., Centrella, J. M., & McMillan, S. L. W., 
1996, Phys. Rev. D, 54, 7261 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v4.0. 



12 



